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We study the spin stiffness of a one-dimensional quantum antiferromagnet in the 
whole range of system sizes L and temperatures T. We show that for integer and half- 
odd integer spin case the stiffness differs fundamentally in its L and T dependence, 
and that in the latter case the stiffness exhibits a striking dependence on the parity 
of the number of sites. Integer spin chains are treated in terms of the non-linear 
sigma model, while half-odd integer spin chains are discussed in a renormalization 
group approach leading to a Luttinger liquid with Aharonov-Bohm type boundary 
conditions. 
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Quantum one-dimensional antiferromagnets have been the subject of intensive studies 
since Haldane |I] conjectured that the spectrum of an integer spin S chain has a finite gap 
even in the absence of any anisotropy, while half-odd integer S chains are gapless. In both 
cases the Neel long-range order of the ground state is destroyed by quantum fluctuations. 
However, the " degree of destruction " is different: for integer S the correlation length is 
finite, which means that the elementary excitations have a gap, while for half-odd integer 
S the correlation length is infinite and excitations are gapless. By now, the presence of 
the Haldane gap for integer S chains is well understood theoretically || and has been 
confirmed in experimental and numerical Q studies. However, all these investigations 
were concentrated on the Haldane gap, i.e., on the energy spectrum itself. Thus a broader 
understanding of the Haldane conjecture is desirable, and the question naturally arises 
whether there are alternative manifestations of the fundamental difference between integer 
and half-odd integer spin chains. Indeed, it is the purpose of this work to provide an 
affirmative answer to this question and to discuss such a particular case in terms of the 
so-called spin stiffness. 

Quite recently, there has been much interest in the spin stiffness (helicity modulus) p s 
of classical Heisenberg ferromagnets 0-0]. p s is defined as a change in the free energy F of 
the magnet when a twist is applied to the spins at the sample boundaries. When thermal 
fluctuations are taken into account, p s is being renormalized with respect to its bare value 
and depends on the scale at which it is probed. Chakravarty has recently shown that 
p s exhibits features which are familiar from the behavior of the electrical conductance of a 
metal in the weak localization regime ||. For instance, in 2d, the mean value of p s depends 
logarithmically on the sample size L, while the rms value of its fluctuations is universal 0. 
This similarity [|Kj makes the spin stiffness an equal member of the family of traditionally 
mesoscopic quantities such as a conductance or a persistent current. Also, the stiffness p s 
serves as a useful (though not perfect) tool for characterizing magnetic long-range order; in 
particular, a vanishing value of p s indicates absence of order [|TT|]. 

Our goal here is to study p s of a quantum one-dimensional antiferromagnet, where fluc- 
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tuations are i) quantum and ii) topologically distinct for integer and half-odd integer S. We 
shall see that the behavior of p s is, indeed, quite different for integer and half-odd integer 
S: p s is renormalized with L in the former case (as it is for a classical 2d ferromagnet), 
whereas it is L-independent (in leading order) in the latter. Moreover, for half-odd integer 
S, p s is shown to exhibit a striking dependence on the parity of the total number N of 
spins. These results should be amenable to a direct check in numerical simulations (see e.g. 
Refs. [f7|Jl2l), and, in particular, could be tested experimentally by measuring the stiffness of 
quasi-one-dimensional antiferromagnets of finite size using similar materials as in Ref. ||. 

We start with the Heisenberg Hamiltonian for a spin chain with nearest-neighbor inter- 
actions 

N 

H = J ex Y,S(n)S(n + l), (1) 

n=l 

where J ex > 0, and we consider the integer S case first. In this case, the long-wavelength 
limit of the partition function Z becomes the (1 + l)d non- linear a model (NLcxM) |T3| with 
Z = J Dri(5(n 2 — 1) exp(— .4.), and the Euclidean action is given by 

A = 7T f dxodxid^n) 2 , p = 0, 1, (2) 
2g Jo Jo 

where g = 2/S, v s = 2SJ ex a is the spin wave velocity (we set k B = h = 1), a is the 
lattice constant, L = clqN, Lt = f3v s is the wavelength of the thermal magnons, and n is the 
slow- varying component of the (staggered) magnetization satisfying the constraint n 2 = 1. 
On the edges of the space-time domain L x Lt the boundary conditions are periodic in xq, 
i.e., n(0, x) = jx(Lt,x), and correspond to a fixed twist of the n-field applied in x-direction, 
i.e., n(x ,0) = (1,0,0) and n(x ,L) = (cos9, sin 9, 0), where 9 is the twist angle. It is 
convenient to use the transformation |]5|,|6|JT^1 n = TZ(9(x))a, where 1Z is the rotation matrix 
about the z-axis by the angle 9(x) = 9x/L, and a satisfies the boundary conditions o\ = 1, 
and a"2 3 = 0, at x = 0, L. The action then takes the form 
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A = Tg J dx[(d,ay + -(1 - ai) 

Q 

+2—(a 1 d x a 2 - a 2 d x a 1 )]. (3) 



We define the spin stiffness p s in units of the velocity 

1 r d2F l m 

where F = —TlnZ ||15|| . With this definition, the bare value of p s in the tree (classical) 
approximation is p° = v s /2g. Corrections due to quantum and thermal fluctuations can be 
found in a loop expansion in g. In one-loop order, only quadratic terms in the action (|3D are 
to be retained, and the third term in Eq. (|3]) reduces to a total derivative d x a2, and, thus, 
vanishes as we shall restrict our consideration to the topological sector with zero winding 
number (Pontryagin index). Performing the functional integration, we get 

where q = (2im/Lj', nm/L), n = 0, ±1, m = 1,2. ..The sum in Eq. (H) can be 
evaluated in three limiting ) L <C Lt (quantum region); b) ao <C Lt <C L (classical 

renormalized region); and c) Lt ao (classical region) . We have 

ln(f 9m /L) /47r, for a) 
^ = \ (e(T) - L)/12Lr, for b) (6) 
(&,-L)/12L T) fore) 

Here, £, qm = aaoexp(7r5') is the correlation length in the quantum region. £(T) = 
3LTln(7^g m /L^)/7r is the classical correlation length ^ = QLx/g renormalized by quan- 
tum fluctuations, a and 7 are (cut-off dependent) non-universal constants of order one. We 
note that case c) agrees with the Id classical NLcrM 0. 

In all the regions, p s goes to zero as the system size L < £ approaches the correlation 
length of the corresponding region ]TB] . This zero value of p s is what a macroscopic system 
is expected to have in the absence of spontaneously broken symmetry []TTj] (the point p s = 
signals the breakdown of the one-loop order expansion). 

The rms fluctuation 5p 2 s = (TLd/d6 2 ) 2 \nZ\g = o is given by 

M = ^_vi- = l®£' fora) m 

^ Ufe) 1 , for b) and c), 



where ((x) is the Riemann ^-function. Contrary to the case of a classical ferromagnet ||, 
the fluctuations are non-universal: they depend both on L and T. Moreover, in the classical 
and classical renormalized regions, the fluctuations are abnormally large (8p s > p s ), and, 
thus, the spin stiffness is not a self-averaging quantity. Finally, we note that the analogues 
of the regions we consider in the quantum NLcrM can also be obtained in the classical 2d 
model, if one considers a rectangular instead of a square system. 

We now turn to the half-odd integer spin case. The effective field-theoretical description 
of the long-wavelength excitations is not of much use here, since the partition function con- 
tains a O-term and contributions from all topological sectors with different winding numbers 
|I|j2]Jl^|, which makes the model hardly tractable. It is believed that the exactly solvable 
case of the spin 1/2 chain reproduces the generic features of all half-odd integer S chains 
|I|j2]jr^|, and we shall consider this case only. By using the Jordan- Wigner transformation 



iKn) = (-l)"exp fiTrXX&O') + S-{n), (8) 

where S± = S x ± iS y , the Hamiltonian ([]]) is mapped on to a system of spinless fermions on 
the lattice, 

N i 

H = J ex £(-- W{n)il>(n + 1) + H. c.} 

n=l 

+ :p{n) ::p(n+l) :), (9) 

where : p(n) := ip^ (ri)i}j{n) — 1/2 is the (Fermi-ordered) density operator. We now have to 
specify the boundary condition for the fermionic operators ip. The quantum generalization 
of the classical boundary conditions for the spin field n, used in the NLcrM treatment of the 
integer S case, is S±(N + 1) = e ±td S±(l), and S Z (N + 1) = S z (l). The boundary condition 
for ip then follows from Eq. (|8|), 

iP{N+l) = e i[7T{NF+N) - e] 4>(l), (10) 

where the number of fermions is Np = N/2, if N is even, and Np = (N + l)/2, if N is 
odd. We have also used the fact that J2S z (n) = 0(1/2) for even (odd) N. The problem 
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defined by Eqs. @ and fllOD is similar to that of spinless electrons on a ring threaded by 
an Aharonov-Bohm flux 9, with the difference that here the boundary conditions depend on 
the parity of N. This parity dependence will result in a striking difference in the behavior 
of p s for even and odd N. 

Finite-size systems of interacting fermions with twisted boundary conditions have re- 
cently been studied in the framework of the Luttinger liquid approach The par- 



ity dependence of the boundary conditions, however, requires a re-examination of the 
bosonization scheme, which we now address. The left- and right-movers are introduced 
by ip(n) = e mkF ip + (n) + e~ mkF ip-(n), where we choose hp = tt/2 for N even and odd. The 
boundary conditions for ip a , where a = ±, take the form 

{e~ l<i ib„(l), for N even 
(11) 
-me^Va(l), for N odd. 

Bosonic fields are introduced by ip a = (27rao) -1 ^ 2 e lv ^^ a , where a = a0 — and dxfl is 



the conjugate momentum of 0. The zero modes of and $ can be chosen in the form [IS 



0o = 0j/v^ + Mv^Fx/L, 
tfo = $m/V^ + (J " 0/n)^(x + \L)/L, (12) 

where J and M are the operators of the topological current and the number of particles 
above the ground state respectively, which satisfy [(f) j, J] = [$m, M] = i. Next, using 



the Baker-Hausdorff formula, we write 

where $ contains contributions from the nonzero modes and from 0j and $m, and is not 
parity dependent; J(M) stands for the eigenvalue of J(M). It is convenient to introduce 
the topological indices k,j and Km, such that nj = 0(1), if J is even (odd); for even N, 
Km = 0(1), if M is even (odd), and, for odd N, k m = 0(1), if M + 1/2 is even (odd). By 
using Eqs. (|TID and (0), we see that kj 5 m must satisfy the following constraints: Kj = 1, 

6 



km = (and vice versa), if N is even, and kj = km, if N is odd. By comparing these 
constraints with the analogous constraints of the fermionic problem [T5J, we can now say 
that the ground state of our spin system is paramagnetic for N even and diamagnetic for iV 
odd. 

The rest of the bosonization procedure is identical to that of Ref. [113, and the bosonized 
Euclidean action takes the sine-Gordon form 

A = C° Jo d 2 x{Md^) 2 + {i/L)^BM 

-(( 7 o/ag)cos(4 v ^0)}, (14) 

where vq = t> s (l+4/7r) 1 / 2 , K = vo/2v s , go = 1/8tt 2 K , and 8 = k j — 9/tt. The bosonic fields 
have been decompactified in the course of the functional integral derivation, and <fi obeys now 
the boundary conditions 4>(xo + k j3vo, x + k±L) = 4>(xq, x) + k y/iin + kiy/7r(2m + k m ), where 
n and m are the winding numbers in xq- and x-directions, respectively. The measure T>cf) of 
the functional integral Z — J T>(f> exp(— .4^) includes the sums over the winding numbers n, m 
and over the topological indices Kj^m- The last term in Eq. (|l~4"D corresponds to Umklapp 
scattering processes between fermions. Since go/K^ is small (~ 0.02), this Umklapp term 
can be treated perturbatively in a standard renormalization group (RG) approach leading 
to the following flow equations ||20|| : 



dg 
dl 



2{K-l)g, g(0)=g 



f = 2vrV, K(0) = K , (15) 

where / = ln(£/ao) with £ = min{L,L-r}. Since we started with the isotropic Heisenberg 
model (JJ), the scaling dimension of the Umklapp term is equal to the critical dimension 
of the model (= 2), i.e., this term is marginally relevant. In this case, the flow proceeds 



along the separatrix between massive and massless phase. On this line the solutions to 
Eqs. flTjl ) are with C » L 



where K* = K(oo) is the fixed point value, and the (non-universal) cut-off L depends 
on (K ,g ) and is of the order a . At the fixed point (£ — > oo), g — and the action 



(|H|) renormalizes to that of a Luttinger liquid with a topological term (oc <9o0), and with 
parameters renormalized through interactions: (Kq,Vo) — > (K*,v*). By comparing with the 
exact Bethe-ansatz solution P^j] , one gets = 1 and t>* = ttv s /2. 

We can now calculate the fixed-point value of the spin stiffness, p*, and its finite size and 
finite temperature corrections. Upon integrating out the zero modes, the twist- dependent 
part of Z becomes 

Z e = E e-^%( Zj ,e- a )e 3 (z M ,e- 4b ), (17) 

KJ,KM 

where 3 (z,q) = J2q n2 e 2mz is the Jacobi 9 function, zj = 7i9 /2, a = -kK*L/L* t , b = 
7tK*Lt/L, zm = 2iKMb, and L* T = (3v* . The results for the spin stiffness take simple forms 
in the limiting cases of low and high temperatures. For L L?, we obtain from (^j) and 
(0) 

-(v*/8K* 2 )(L*/L), for iV even 
v / ; ' (18) 

v*/AttK*, for odd, 
while for L ^> L* T we get 

= (-l) N+1 2(L/L* T )e- nxL/L T, (19) 

where x — K* + l/4i^*. To obtain the value of p s away from the fixed point, we go back 
to the full action QUI) and replace (K ,g ) by (K,g) from Eq. (|T6|) , treating the deviations 
from the fixed-point values as perturbations pT|-|23H . In first order, the Umklapp term gives 
no contribution, while the perturbation in K — K* leads to: 

, ( 1 + l/KML/Lo)], for L « L* T 
p 8 = p*x\ (20) 
[ exp (3Ltt/[8L^ \n(T /T)}) , for L > Ly, 

where a at = 1 for even N, and a at = 2 for odd iV; the cut-off T is of the order of 
v*/Lq. The last equation is valid for T < T* < T , where L^(T*) = ao- The L- and 
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T-dependent corrections to p* resulting from the perturbation of the fixed-point action 
by the marginally irrelevant operator are larger than those coming from the expansion of 
Eq. ( [TTD in (L/L^)^ at the fixed point. In particular, the exponential dependence of p* 
on K* in the high-temperature regime results in a significant T-dependent renormalization. 
This renormalization may be conjectured to remain significant in the intermediate regime 
L ~ L* T as well. We also note that the Umklapp processes lead to the breakdown of the 
single-parameter scaling of p s : the latter scales with L/Lj, at the fixed-point but acquires 
additional L/Lq and T/Tq scalings away from the fixed point (for a rough estimate, Lq ~ do, 
and T ~ J ex)- This breakdown can be detected in numerical and real experiments. The 
fluctuations in p s can be calculated along the same lines as for p s itself, and, in marked 
contrast to the integer S case, turn out to be exponentially small for all L/L?- 

The negative value of p s for even N simply reflects the fact that in this case the free 
energy has a maximum at 9 = 0, and an arbitrarily small twist drives the system out of this 
state. Analysing p s at finite 9 (and low temperatures) we can see that p s vanishes at some 
9* ~ LjL* T <c 1 and then remains positive for all 9, thus exhibiting a crossover from the 
paramagnetic to the diamagnetic regime. The parity effects in the spin stiffness are quite 
similar to that in persistent currents of electronic systems P5|pB| , |rS|| ; in particular, the result 
obtained above can be checked without approximations for the special case of the XF-model 
by mapping it on to the exactly solvable problem of free fermions p7j . 

This work was supported by the NSERC of Canada (D.L. and D.L.M.) and in part by 
the US NSF through Grant No. DMR-89-20538 (D.L.M.). 
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